Orientation-Dependent Transparency of Metallic Interfaces 
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As devices are reduced in size, interfaces start to dominate electrical transport making it essential 
to be able to describe reliably how they transmit and reflect electrons. For a number of nearly 
perfectly lattice-matched materials, we calculate from first-principles the dependence of the interface 
transparency on the crystal orientation. Quite remarkably, the largest anisotropy is predicted for 
interfaces between the prototype free-electron materials silver and aluminium for which a massive 
factor of two difference between (111) and (001) interfaces is found. 

PACS numbers: 71.15-m, 72.10.Bg, 72.25.Ba, 73.40.Cg, 73.40-c, 75.47.De 



A recurring theme in condensed matter physics in the 
last twenty years has been the discovery of new physical 
effects and properties in systems with reduced dimen- 
sions; the prospect of exploiting these effects and prop- 
erties in logical processing, sensing and storage devices 
is an important driving force behind nano-science and - 
technology. In semiconductors, the electronic structures 
of the electrons responsible for conduction can be de- 
scribed using simple models. The same is not true of the 
ferromagnetic transition metals which form the basis for 
magnetoelectronics. It is the non-trivial spin-dependence 
of the transmission and reflection of electrons at magnetic 
interfaces which provides the key to understanding phe- 
nomena such as oscillatory exchange coupling, giant- and 
tunneling- magnetoresistance, spin transfer torque, spin- 
pumping and spin injection F° r well-studied material 
combinations such as Co/Cu and Fe/Cr, modest spin- 
dependence of the interface transmission 0, |j| Q of the 
order of 10-20% is sufficient to account for experimental 
observations 0- 

However, the confrontation of theory and experiment 
just referred to is at best indirect and model-dependent. 
Even though the theory of transport in small structures 
is formulated in terms of transmission and reflection ma- 
trices , measuring interface transparencies directly has 
proven quite difficult pj. To identify interfaces suitable 
for experimental study, we have undertaken a systematic 
materials-specific study of the orientation dependence of 
the interface transmission between pairs of isostructural 
metals whose lattice constants match within a percent or 
so in the hope that it will prove possible to grow such 
interfaces epitaxially. 

One of the metal pairs we studied was Al/Ag. Both 
metals have the fee crystal structure and their lattice con- 
stants arc matched within 1%. Aluminium is a textbook 
8] example of a system well described by the (nearly) free 
electron model. Silver, also usually assumed to be a free 
electron-like material, is a noble metal with high conduc- 



tivity which is frequently used for electrical contacting. 
We found that in spite of the simplicity of both metals' 
electronic structures the transmission through Al/Ag in- 
terfaces can differ quite significantly from the predictions 
of the free electron model. In particular, between (111) 
and (001) orientations we find a factor 2 difference in in- 
terface transmission for clean Al/Ag interfaces. For free 
electrons the anisotropy should vanish. Our result is in- 
sensitive to interface disorder. We identify a new factor 
responsible for this difference which is not related to the 
standard velocity- 0, El or symmetry-mismatch 0, ll2T ] 
mechanisms. 

A free electron description of interface scattering, in 
which the effect of the crystal potential on transport is 
completely neglected, underlies the BTK theory [1J| used 
to interpret UE3 Andreev reflection (AR) experiments. 
Point contact AR experiments are increasingly used to 
identify the pairing symmetry of superconductors and, 
in the field of magnetoelectronics, to determine the po- 
larization of magnetic materials [liL Ilij . Our finding 
that the electronic structure can have such a large effect 
on the interface transmission, implies that experiments 
should be analysed using more sophisticated models. 

Our study was based upon first-principles calculations 
of the interface electronic structure performed within the 
framework of density functional theory (DFT) and the 
local spin density approximation (LSD A). Bulk and in- 
terface potentials were determined self-consistently us- 
ing the tight binding linearized muffin tin orbital (TB- 
LMTO) [Tg surface Green's function method 0. We 
assumed common lattice constants for both metals of a 
given structure e.g. clai — a-Ag — 4.05 A. The potentials 
obtained in this way were used as input to a TB-MTO 
wavefunction-matching 0, calculation of the trans- 
mission and reflection coefficients between Bloch states 
on either side of the interface. The efficiency of this ap- 
proach is such that interface disorder can be modelled 
using large lateral supercells. For disordered systems, 



the potentials were calculated using the layer CPA (co- 
herent potential approximation) . The results of the cal- 
culations for a number of lattice-matched materials are 
summarized in Table [I] 

The Sharvin conductances, Ga and Gb, reported in 
the third and fourth columns of Tableware proportional 
to the number of states at the Fermi level propagating in 
the transport direction. They are properties of the bulk 
materials which are determined by the area of the Fermi 
surface projections and are a measure of the current- 
carrying capacity of the conductor in the ballistic regime. 
The largest intrinsic orientation dependence, seen to be 
about 13 %, is found for Mo; for Al and Ag, respectively, 
it is less than 8% and 5%. 

The interface transmission in column five of Table [I] 
is expressed as a conductance, Ga/b = e 2 /h T^ v , 
where T^„ is the probability for the incoming state v 
in material A to be transmitted through the interface 
into the out goin g state [i in material B. For most pairs 
of materials |l8l|. the orientation dependence of Ga/b is 
modest (~ 15% for Mo/W) and the interface conduc- 



TABLE I: Sharvin conductances and interface transmissions 
in units of 10 lb Q.~ 1 m~ 2 , interface resistances SR 0, 0] for 
ideal (and, in brackets, for disordered) interfaces in units of 
W~ 15 Q,m 2 . S is the area of the sample for which R is mea- 
sured. Interface disorder was modelled in 10 x 10 lateral su- 
percells with two layers of 50-50 alloy. The largest uncer- 
tainty between different configurations of disorder is about 
2.3%. The values given are for a single spin. For the pairs of 
materials and orientations indicated by a (*), comparison of 
the interface resistances shown in the last two columns with 
experimental values extracted from measurement on multi- 
layers by the MSU group 5, 7) yields reasonable quantitative 
agreement |2l I tll l2|. 



A/B G A G b G a/b 2SR 

Al/Ag (TTT) 0.69 0.45 0.41 (0.36) 0.64 (0.92) 

&fcc = 4.05 A (110) 0.68 0.47 0.30 (0.32) 1.60 (1.39) 

(001) 0.73 0.45 0.22 (0.24) 2.82 (2.37) 

Al/Au (TIT) 0.69 0.44 0.41 (0.35) 0.60 (0.99) 

a /cc = 4.05 A (001) 0.73 0.46 0.24 (0.26) 2.37 (2.14) 

Pd/Pt (111) 0.62 0.71 0.55 (0.54) 0.30 (0.33) 

a/ee = 3.89 A (001) 0.58 0.70 0.52 (0.51) 0.37 (0.39) 

W/Mo (001) 0.45 0.59 0.42 (0.42) 0.42 (0.42) 

abcc = 3.16 A (110) 0.40 0.54 0.37 (0.38) 0.52 (0.47) 

Cu/Co (111)* 0.56 0.47 0.43 (0.43) 0.34 (0.35) 

majority (001) 0.55 0.49 0.46 (0.45) 0.26 (0.27) 

a /ee = 3.61 A (110) 0.59 0.50 0.46 (0.46) 0.35 (0.35) 

Cu/Co (111)* 0.56 1.05 0.36 (0.31) 1.38 (1.82) 

minority (001) 0.55 1.11 0.32 (0.32) 1.79 (1.79) 

a /ee = 3.61 A (110) 0.59 1.04 0.31 (0.35) 1.89 (1.55) 

Cr/Fe (111) 0.61 0.82 0.27 (0.31) 2.22 (1.84) 

majority (001) 0.64 0.82 0.11 (0.25) 7.46 (2.55) 

a, bcc = 2.87 A (110)* 0.59 0.78 0.22 (0.27) 3.04 (2.18) 

Cr/Fe (111) 0.61 0.41 0.34 (0.34) 0.93 (0.95) 

minority (001) 0.64 0.46 0.35 (0.35) 0.98 (0.95) 

& bcc = 2.87 A (110)* 0.59 0.40 0.32 (0.32) 1.03 (1.06) 




FIG. 1: Top row: Fermi surface projections for (a) Ag, (b) 
Al and (c) transmission probabilities in 2DBZ for (111) orien- 
tation. Middle row: Same for (001) orientation. The colour 
bars on the left indicate the number of scattering states in the 
leads for a given two dimensional wavevector ky . The trans- 
mission probabilities indicated by the colour bars on the right 
can exceed 1 for kys for which there is more than one scatter- 
ing state in both Ag and Al. Bottom row: Fermi surfaces of 
(g) Ag and (h) Al, (i) the interface adapted BZ for (001) and 
(111) orientations. The vertical dashed line in (c) and on the 
yellow plane in panel (i) indicate the cross-section used in the 
left-hand panel of Fig. [5] while the vertical dotted line in (f) 
and on the blue plane in panel (i) indicate the cross-section 
used in the right-hand panel. 

tance itself tends to be slightly smaller than the lower 
of the two Sharvin conductances. For these systems the 
behaviour of the transmission appears to be determined 
by the projection of the Fermi surfaces. However, this is 
not so for Al/Ag and Al/Au interfaces. Here we observe 
a large anisotropy in the transport properties. The fac- 
tor 2 difference in transmission between (111) and (001) 
orientations results in a factor 4 difference between 
interface resistances estimated using the method of |2j, |4( . 

The transmission probability for the (111) and (001) 
orientations is plotted in Figs.^ and Of as a function of 
the wave- vector component parallel to the interface, k|| , 
within the 2D interface Brillouin zones (BZ). A qualita- 
tive difference between the two orientations can be ob- 
served. In the (111) case, the transmission is almost uni- 
formly high wherever there are states on both sides of the 
interface. The (001) orientation exhibits more variation 
with high transmission in the central and outer regions of 
the 2D BZ but much lower in a ring-shaped region in be- 
tween. The presence of this "cold ring" is the reason why 
the total transmission is lower for the (001) orientation. 
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FIG. 2: Intersection of a (110) plane with the Al Fermi surface 
and with the interface adapted BZs indicated in Fig.0 (where 
the meaning of the dashed and dotted lines is explained) . The 
labelled dots indicate the positions of the RL sites with coor- 
dinates given in units of 2n/a. The red (blue) lines indicate 
regions of high (low) transmission. 



Explaining the transparency anisotropy of Al/Ag inter- 
faces requires finding an explanation for the low trans- 
mission values in this region of the 2D BZ. 

Two mechanisms are usually taken into account when 
analysing the scattering at perfect interfaces. The first, 
velocity mismatch, is the modulation of the transmis- 
sion by a factor reminiscent of the free electron for- 
mula for the transmission through a potential step: T = 
4vlVr/(vl + vr) 2 where vl/r are the components of 
the Fermi velocities in the transport direction on the 
left and right sides of the interface. This modulation 
is indeed present in our calculated transmissions but its 
effect tends to be noticeable only when one of the ve- 
locities is almost vanishingly small. Naive application of 
the free electron formula yields uniformly good trans- 
mission |20j independent of the orientation. Symme- 
try mismatch, the second mechanism, can suppress the 
transmission between states of incompatible symmetries 
{e.g. even vs. odd etc.). Examination of the eigenvec- 
tors demonstrates that this is not the case for the Al/Ag 
system. For example, states on both sides of the inter- 
face, with k|| along the vertical dotted line in Fig.^, are 
even under reflection in the plane denned by this line and 
the (001) transport direction. Their orbital composition 
(s ,p y ,p z ,d yz ,d,3 Z 2 _ r 2 ,d x 2 _ y 2 where the y axis is parallel to 
the dotted line and z is the transport direction) is essen- 
tially the same for both materials. The same holds for 
states along other symmetry lines/planes and general k|| 
points (in the sense of orbital composition). The origin 
of the "cold ring" must be sought elsewhere. 

In spite of the failure of the free electron transmis- 
sion formula, this simple model serves as a useful starting 
point for analysing the Fermi surface (FS) topologies. In 
the simplest possible approach, we model the FS of Ag 
(shown in Fig. Q^) as a sphere which fits into the first 
BZ. A larger sphere, accommodating three electrons, is 
needed for trivalent Al. In an extended zone scheme, 



conservation of momentum parallel to the interface dic- 
tates that the transmission through a specular interface 
is non-zero only between states with the same values of 
k|| ; these are the k|| -vectors belonging to the region where 
projections of the Fermi spheres on a plane perpendicu- 
lar to the transport direction overlap. For systems with 
lattice periodicity, we must use a downfolded FS, with 
fragments of the original FS sphere back-translated into 
the 1st Brillouin zone, a procedure which can be realized 
geometrically by placing spheres accommodating three 
electrons on reciprocal lattice (RL) sites and then only 
considering the fragments in the first BZ. Examination 
of the FS of Al calculated from first-principles (Fig. HJi) 
and its cross-section (Fig. EJ) reveals that, in spite of its 
apparent complexity, it remains essentially (piecewise) 
spherical. For some values of k|| (see Figs, and^), 
Al can now have more than one propagating state. Nev- 
ertheless, in the free electron limit, the downfolded states 
are strictly orthogonal to the states in Ag and the total 
transmission is unchanged. For a reduced zone scheme, 
we formulate the following rule: The transmission be- 
tween states in two (nearly) free electron materials which 
have the same kn , but originate from reciprocal lattice 
sites whose parallel components do not coincide, vanishes 
in the free electron limit and is expected to be strongly 
suppressed for nearly free electron like materials. 

Obviously, the truly free electron system can not ex- 
hibit anisotropy. However, in the presence of the periodic 
potential the original, piecewise-spherical Fermi surface 
and consequently the transmission is going to be mod- 
ified. Firstly, since the wave functions are no longer 
pure plane waves, the strict orthogonality of the down- 
folded states is relaxed and the transmission can assume 
finite although typically small values (hence suppressed 
instead of zero in the above rule). Secondly, the shape 
of the Fermi surface changes with the modifications be- 
ing strongest in the vicinity of RL planes where, for Al, 
we observe the opening of gaps between previously con- 
nected fragments. The anisotropy is mostly related to 
this second effect. 

In Fig. |2 we show the intersection of the Al FS with a 
(110) plane. The two plots are rotated so that the vertical 
axis in Fig. is the [111] direction while in Fig. [2t> it is 
[001]. In both cases the positions of the nearest RL sites 
(on which spheres are centred) are shown together with 
the cross-section through the relevant interface-adapted 
Brillouin zone, which is different for each orientation; 
see Fig. We can now readily identify spheres from 
which various fragments of the Fermi surface originate 
and mark those fragments with positive (upward) veloc- 
ities, according to the rule given above, as having high 
(red) or low (blue) transmissions. In the (001) case, the 
"high" fragments originate from (0,0,0) and (0,0,-2) cen- 
tred spheres. Comparing Figs. ^ and 03, we note that 
the position of the gaps opened between these spheres by 
Bragg reflection on the (001) and (001) planes coincides, 
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in projection along the [001] direction, with the position 
of the "cold ring" in Fig. The other states present in 
this region originate from (1,-1,-1) (and equivalent) cen- 
tred spheres, are therefore nearly orthogonal to states in 
Ag centered on (0,0,0) and so have low transmission. In 
the (111) case however, the large fragments of FS belong- 
ing to the same (1,-1,-1) sphere have high transmissions 
(Fig. Et) an d dominate transport. In addition, the effect 
of gap-opening is reduced in this orientation because of 
the rotation. Combination of these two factors results in 
the almost uniformly high transmission seen in Fig. 

We can now finally identify the origin of the trans- 
mission anisotropy for Al/Ag interface. It stems from 
two factors: (i) the near orthogonality of the downfolded 
Al states to those belonging to the simple Ag sphere and 
(ii) the gaps opened in the continuous free electron Fermi 
surface by the periodic potential. The latter factor is of 
course related to the symmetry of the underlying crys- 
tal lattice and directly responsible for the introduction 
of the orientation dependence. For Al/Au interfaces, the 
interface transmissions and resistances are very similar 
to the Al/Ag case. 

The orientation dependence of the interface transmis- 
sion of six metal pairs with the same structure and lattice 
parameter was calculated. For fee Ag/Al a factor two 
difference between the (111) and (001) orientations was 
found and explained within the free electron model. The 
predicted anisotropic interface resistance and Andreev re- 
flection (not shown) are not very sensitive to interface 
disorder and should be observable experimentally. 
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